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We present a novel comprehensive first-principles theoretical study of the electronic properties and 
relative stabilities of edge-oxidized zigzag graphene nanoribbons. The oxidation schemes considered 
include hydroxyl, carboxyl, ether, and ketone groups. Using screened exchange density functional 
theory, we show that these oxidized ribbons are more stable than hydrogen-terminated nanoribbons 
except for the case of the etheric groups. The stable oxidized configurations maintain a spin-polarized 
ground state with antiferromagnetic ordering localized at the edges, similar to the fully hydrogenated 
counterparts. More important, edge oxidation is found to lower the onset electric field required 
to induce half-metallic behavior and extend the overall field range at which the systems remain 
half-metallic. Once the half-metallic state is reached, further increase of the external electric field 
intensity produces a rapid decrease in the spin magnetization up to a point where the magnetization 
is quenched completely. Finally, we find that oxygen containing edge groups have a minor effect on 
the energy difference between the antiferromagnetic ground state and the above-lying ferromagnetic 
state. 



Low-dimensional carbon structures such as fullerenes 1 
and carbon nanotubes 2 (CNTs) are promising candidates 
for building blocks of future nanoelectronic and nanome- 
chanical devices. 3 ^ Made of a unique hexagonal carbon 
lattice confined to a quasi one-dimensional (ID) tubular 
structure, CNTs may be either semiconducting or metal- 
lic, depending on their diameter and chirality. 5 Combined 
with their ballist io&LSa2 n 10 1 1 1 electronic transport charac- 
teristics, this opens exciting possibilities for the design 
of novel electronic components and interconnects with 
nanometer scale dimensions. 

Recently, a new type of graphene-based material was 
experimentally realized. 12 Shaped as narrow stripes cut 
out of a single (or a few) layer(s) of graphite, these 
elongated materials were named graphene nanoribbons 
(GNR). Since they share the same hexagonal carbon lat- 
tice structure, GNRs and CNTs exhibit many similarities 
with respect to their electronic properties. However, the 
planar geometry of the ribbons allows for the applica- 
tion of standard lithographic techniques for the flexible 
design of a variety of experimental devices i 12 i 13 i 14 i 15 i 16 
The utilization of these established fabrication methods 
suggest the possibility of a controllable and reproducible 
fabrication of carbon based electronic components at the 
nanometer scale. 

Despite the aforementioned similarities, there is a dis- 
tinct difference between CNTs and GNRs. Unlike CNTs, 
GNRs present long and reactive edges prone to localiza- 
tion of electronic edge states and covalent attachment of 
chemical groups^ that can significantly influence their 
electronic properties. The importance of edge states 
was demonstrated for the case of arm-chair CNTs which 
present a metallic and non-magnetio 3 -^ character in their 
tubular form, but when unrolled into the corresponding 
zig-zag GNRs, they were predicted to become semicon- 
ducting^^ with a spin polarized^^^ 2 ^ 3 .* 2 ^ 5 . ground 
state. This ground state is characterized by opposite spin 
orientations of localized electronic states at the two edges 
of the GNR, which couple through the graphene back- 
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FIG. 1: Optimized unit cell geometries of different edge oxi- 
dation schemes studied in this work. Color code: red - oxygen 
atoms, gray - carbon atoms, white - hydrogen atoms. 



bone via an antiferromagnetic (AF) arrangement of spins 
on adjacent atomic sites. 

Recently, Son et al 19 have shown that upon the appli- 
cation of an electric field, an opposite local gating effect 
of the spin states on the two edges of the ribbon may oc- 
cur. The in-plane field (perpendicular to the main axis 
of the ribbon) drives the system into a half-metallic state 
where one spin flavor exhibits a metallic behavior, while 
the opposite spins experience an increase in the energy 
gap. Apart from the interesting interplay between the 
electric field and the spin degree of freedom, this phe- 
nomenon is important from a technological standpoint 
since such a system could serve as a perfect spin filter in 
future nanospintronic devices. 

Most of the theoretical studies on 

QNRg 18.20.21.22.23.24.25.26.27.28.29.30.31.32.33.34.35.36 

eluding those on the half- metallic systems , 19 i 37 address 
either bare or hydrogen terminated ribbons. Neverthe- 
less, during standard GNRs fabrication processes, it 



2 



is commonly assumed that the ribbons edges become 
oxidized j 38 i 39 i 40 i 41 Since many of the unique properties 
of GNRs are associated with edge states, edge chemistry 
in general and oxidation in particular, can significantly 
alter the electronic properties of the ribbons. It is the 
purpose of this letter to study, for the first time, the 
influence of edge oxidation on the relative stability, the 
electronic properties, and the half-metallic nature of 
zigzag graphene nanoribbons, from first- principles. 

To this end, we study six different oxidation schemes 
of a 1.8 nm wide zigzag GNRs including hydroxylation, 
carboxylation, ketonation, and etheration, as shown in 
Figure [TJ All the calculations presented in this work 
were carried out using periodic boundary conditions as 
implemented in the development version of the Gaus- 
sian suite of programs i 4 ^ 4 - 3 Fully relaxed geometries were 
obtained using the PBE realization of the generalized 
gradient approximation 44,45 and the polarized 6-3 1G** 
Gaussian basis set 46 for the non-magnetic (closed-shell) 
state. The electronic structures were then re-evaluated 
using the screened exchange hybrid density functional, 
HSE, 47,48 which has been tested in a wide variety of 
materials and has been shown to accurately reproduce 
experimental band gap a 49 i 5Q and first and second op- 
tical excitation energies in metallic and semiconduct- 
ing SWNTs j 51 i 52 Furthermore, the inclusion of short- 
range exact-exchange in the HSE functional makes it 
suitable to treat electronic localization effect o 53 i 54 i 55 i 56 i 57 
which are known to be important in this type of materi- 
als . 18 ' 19 ' 2Q ' 21 ' 22 ' 23 ' 24 ' 25 - 26 ' 27 ' 28 ' 37 - 58 ' 59 ' 6Q ' 61 ' 62 - 63 

We start by studying the relative stability of the dif- 
ferent oxidized ribbons. As these structures have differ- 
ent chemical compositions, the cohesive energy per atom 
does not provide a suitable measure for the comparison of 
their relative stability. Therefore, we adopt the approach 
customary used in binary phase thermodynamics to ac- 
count for chemical composition and utilized previously to 
qualitatively analyze the relative stability of endohedral 
silicon nanowires and arm- chair GNRs. 31 Within this 
approach one defines a Gibbs free energy of formation SG 
for a GNR as: 

SG(x) = E(x) - XhI^h ~ Xol^o ~ XcVc, (1) 

where E(x) is the cohesive energy per atom of the GNR, 
Xi is the molar fraction of atom i (i=C,0,H) in the rib- 
bon, satisfying the relation T>iXi = 1, and \ii is the chem- 
ical potential of the constituent i at a given state. We 
choose ji H as the binding cohesive energy per atom of 
the singlet ground state of the H2 molecule, jjjQ as the 
cohesive energy per atom of the triplet ground state of 
the O2 molecule, and [i c as the cohesive energy per atom 
of a single graphene sheet. This definition allows for a 
direct energy comparison between oxidized nanoribbons 
with different compositions, where negative values rep- 
resent stable structures with respect to the constituents. 
It should be stressed that this treatment gives a qualita- 
tive assessment of the relative stability while neglecting 
thermal and substrate effects and zero point energy cor- 
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FIG. 2: Ground state relative stabilities of the different ox- 
idation schemes studied (see Fig. [1]) obtained via Eq. [1] at 
the PBE and HSE/6-31G** levels of theory. Negative values 
indicate stable structures with respect to the constituents. 



rections. 

In Figure [2] we present the relative stability of the 
ground states of the different oxidation schemes stud- 
ied. The fully hydrogenated ribbons, as well as the two 
etherated ribbon configurations, are found to be less sta- 
ble than their corresponding constituents. This is con- 
sistent with previous calculations on bare and hydrogen 
terminated GNRs. 31 Nevertheless, the chemical passiva- 
tion via hydroxylation, ketonization and carboxylation 
leads to considerable energetic stabilization of the struc- 
ture of the ribbons. We find that the most stable struc- 
ture corresponds to the fully hydroxylated ribbon. This 
enhanced stability is attributed to the hydrogen bonds 
formed between adjacent hydroxy 1 groups. 

We now turn to study the electronic properties of the 
four most stable structures. In Figure [3] we show the 
ground state spin density of the fully hydrogenated, ke- 
tonated, hydroxylated and carboxylated structures, ob- 
tained with the HSE functional^ Similar to the fully 
hydrogenated GNR 19 , all the depicted oxidized struc- 
tures exhibit a spin polarized ground state where the spin 
magnetization on the opposite edges of the ribbons are 
aligned anti-parallel. For the hydroxylated and the car- 
boxylated ribbons, the presence of the oxidation groups 
slightly changes the absolute value of the spin polariza- 
tion at the edge carbon atoms. However, the overall 
magnetization density resembles that of the fully hydro- 
genated GNR. A different picture arises for the ketonated 
system, where the oxygen p-electrons participate in the 
7r system of the ribbon and therefore present consider- 
able spin polarization, which results in a qualitatively 
different spin density map compared to that of the fully 
hydrogenated system. 

In Table [J we present the bandgaps and the energy 
differences per unit cell between the AF ground state 
and the above-lying ferromagnetic (FM) state (where the 
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FIG. 3: Map of the HSE ground state spin densities 0.16 A 
above the surface of the ribbon for the four most stable ox- 
idation schemes studied. The fully hydrogenated system is 
also included for comparison. Color code: Red a spin den- 
sity; blue (3 spin density. A scheme of the different oxidation 
groups is presented for clarity. 



FIG. 4: Electric field effect on the spin polarized bandgap 
of oxidized GNRs. Red and blue triangles stand for the two 
different spin polarization bandgaps. The values are calcu- 
lated using the HSE functional and the 3-21G basis set.— A 
scheme of the direction of the applied electric field is presented 
for clarity, where X and Y represent the oxidizing groups. 



TABLE I: HSE/6-31G** AF ground state bandgaps (eV) and 
energy differences between the FM and AF states (meV/unit 
cell) of the studied GNRs. 



Eg E F m — Eaf 
Hydrogenated 1.05 42 
Hydroxylated (I) 0.99 39 
Hydroxylated (II) 0.90 38 
Carboxylated 0.77 33 
Ketonated 0.03 5 



edges of the ribbon bare parallel spin polarization) for the 
hydrogenated, and stable oxidized structures. All stud- 
ied oxidized structures, except for the ketonated ribbons, 
have a bandgap comparable to that of the fully hydro- 
genated ribbon, indicating that their electronic character 
is only slightly changed upon oxidation. Furthermore, 
the energy difference between the AF ground state and 
the FM state is above room temperature and similar to 
that of the hydrogenated ribbon. The ketonated ribbon, 
on the other hand, exhibits a vanishing bandgap and a 
very small energy difference between the AF ground state 
and the FM state, suggesting that it will be extremely 
difficult to observe the half-metallic nature of the keto- 
nated ribbons even at low temperatures. 

The analysis presented above suggests that, with the 
exception of the ketonated systems, all oxidation schemes 
have little effect on the electronic character of the GNRs 



and therefore, one would expect these systems to behave 
as half-metals under the influence of an external elec- 
tric field. To verify this assumption we have calculated 
the bandgap of the a and (3 spin channels as a function 
of the intensity of a transverse in-plane, uniform, and 
static electric field. In Figure [H we plot the spin resolved 
bandgaps as a function of the field intensity for the fully 
hydrogenated, hydroxylated, and carboxylated ribbons. 
In the absence of an electric field, the a and (3 bandgaps 
are degenerate for all the studied systems. Upon the 
application of a low intensity electric field, a splitting oc- 
curs where the bandgap of one spin flavor increases and 
that of the opposite spin flavor decreases. 19 The bandgap 
splitting increases monotonically with the field intensity 
up to a point where the system becomes half- metallic. 
Further increase in the external electric field intensity re- 
sults in a decrease of the bandgap splitting up to a point 
where the systems become non- magnetic. 

It is interesting to note that when the edges of the rib- 
bon are fully or partially hydrogenated, the field intensity 
needed to switch the system to the half-metallic regime 
is 0.4 V/A and the range at which the half-metallic be- 
havior is maintained is 0.3 V/A. Nevertheless, when the 
edges are fully oxidized, the system turns half-metallic at 
a lower field intensity (0.2 V/A) and the range of half- 
metallic behavior is doubled to 0.6 V/A. As shown pre- 
viously in Ref. [l9|, the intensity of the field needed to 
achieve half-metallicity decreases with the ribbon width. 
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FIG. 5: Largest magnetization splitting between the two 
edges of the ribbon as a function of the electric field intensity. 
Calculations performed on the HSE/3-21G level of theory. 



For a given width, the decrease of the field intensity 
needed to obtain the half- metallic behavior combined 
with the enhancement of the range at which this behav- 
ior is maintained, implies that edge oxidation is critical 
for the fabrication of robust and chemically stable spin 
filter devices out of zigzag GNRs. 

It is relevant to investigate how the spin magnetization 
depends on the intensity of the applied electric field. To 
this end, in Figure [5] we present the maximum magne- 
tization splitting as a function of the electric field, i.e. 
the difference between the maximum a and f3 Mulliken 
atomic spin densities. 

For the four systems studied, at the low electric field 
regime, the field has minor influence on the magnetiza- 
tion splitting. Once the systems become half-metallic, 



the magnetization drops sharply with the applied field 
up to a point where it fully quenches. An explanation 
for this behavior can be found by realizing that once the 
systems become half- metallic, the electric polarizability 
of the metallic spin channel increases considerably, allow- 
ing for a more effective spin compensation between the 
opposite edges of the ribbon. This results in an overall 
quenching of the spin magnetization. The electric field 
above which all the systems become diamagnetic is 0.7- 
0.8 V/A for ribbons of about 1.8 nm wide. 

In summary, we have presented a detailed study of 
the electronic properties and relative stabilities of sev- 
eral edge-oxidized zigzag GNRs. The structure of the 
oxidized ribbons is found to be stabilized with respect to 
the fully-hydrogenated counterparts except for the case of 
the ether ic groups. Notably, all the stable oxidized struc- 
tures studied maintain a spin polarized ground state with 
antiferromagnetic ordering. The above-lying ferromag- 
netic arrangement is expected to become accessible only 
beyond room temperature. Apart from the ketonated 
ribbons, all oxidized systems present a lower onset field 
and a wider electric field range for which half-metallic 
behavior is predicted. This suggest that edge oxidation 
is important for the design of efficient and robust spin- 
tronic devices based on GNRs. We have also shown that 
once the half-metallic state is reached, further increase of 
the external electric field produces a rapid decrease in the 
spin magnetization up to a point where all the systems 
become non-magnetic. 
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